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Abstract. Ishimori equation is a (2 + 1) dimensional generalization of the (1 + 1) 
dimensional integrable classical continuous Heisenberg ferromagnetic spin equation. 
The richness of the coherent structures admitted by Ishimori equation I such 
as dromion, lump and rationally- exponentially localized solutions, have been 
demonstrated in the literature through d technique and binary Darboux transformation 
method. To our knowledge Hirota's method had been adopted to construct only 
the vortex solutions of Ishimori equation II. For the first time, the various types of 
localized coherent structures mentioned above have been constructed in this paper for 
the Ishimori equation I using the Hirota's direct method. In particular we have brought 
out the significance of boundaries and arbitrary functions in generating all these types 
of localized structures and proved that the absence of such boundaries leads only to 
line soliton solutions. 
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1. Introduction 

A (2 + 1) dimensional integrable generalization of the (1+1) dimensional integrable 
Heisenberg ferromagnetic spin equation (isotropic Landau-Lifshitz equation) S t (x,t) = 
SAS XX [1] was introduced by Ishimori in 1984 [2] to explain the dynamics of the classical 
spin system on a plane. Its form is 

S t (x, y, t) = S A (S xx + <y 2 S yy ) + <p y S x + (f> x S y , (1.1a) 

4>XX - 0- 2 4>yy = ~2a 2 S- S x A Sy, (1.16) 

where S = (Si, S2, S3) is the three dimensional spin unit vector (S 2 = 1), <p(x, y, t) is a 
scalar field and a 2 = ±1. If a 2 — +1, eq. (1.1) is referred to as Ishimori equation I (IE 
I) and if a 2 — —1, it is referred to as Ishimori equation II (IE II). An important feature 
associated with eq. (1.1) is the existence of nontrivial topological invariant known as 
topological charge defined as 

Q = ^J J ' S-S* A Sy dxd V ( L2 ) 

and the solutions of eq. (1.1) are classified in terms of the integer values of Q. Just 
like the (1+1) dimensional integrable spin system is geometrically equivalent [1] to the 
nonlinear Schrodinger equation through a moving space curve formalism, eq. (1.1) is 
geometrically equivalent to the Davey-Stewartson equation through a moving surface 
formalism [3]. 

The initial value problems of both IE I and IE II have been analysed by the d 
and nonlocal Riemann-Hilbert problem methods, respectively, in [4] by Konopelchenko 
and Matkarimov. For stationary boundaries, the initial boundary value problem 
for the IE I has been studied in [5] and three different types of localized solutions 
(soliton-soliton(ss), soliton-breather(sb), breather-breather (bb)) have been presented. 
The line solitons admitted by both IE I and IE II have been presented in [6]. The 
localized coherent structures for the IE I have been analysed in [7] for time dependent 
boundaries and the solutions such as rationally localized soliton, exponentially localized 
soliton and rationally-exponentially localized soliton have been reported using inverse 
scattering transform (1ST) method. IE I has also been analysed through binary Darboux 
transformation method and different types of solutions have been constructed in terms 
of grammian determinants in [8]. Curiously, the Hirota's bilinearization method, which 
is one of the celebrated direct methods and applicable to almost all integrable soliton 
equations in (1+1) and (2+1) dimensions, has not been applied so far (as far as our 
knowledge goes) to obtain localized solutions for IE I, though it has been used by 
Ishimori himself [2] to obtain vortex solutions to IE II. The difficulty probably lies 
in introducing the boundaries in the bilinearized form appropriately for IE I. In this 
paper, we have successfully obtained all types of localized structures of IE I reported 
in [5-7] through the Hirota's bilinearization technique and presented the different types 
of localized structures admitted by it for the case of both time dependent and stationary 
boundaries. 
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The paper is organised as follows: In section 2, eq. (1.1) is bilinearized in laboratory 
coordinates through stereographic projection and Painleve property and multiline 
soliton solutions for both IE I and IE II are presented for completeness. In section 
3, the bilinearized version of IE I in terms of light cone coordinates is presented and the 
general form of the solution is given. In particular, we point out how boundaries can 
be introduced explicitly into the bilinearized version of eq. (1.1). The role of linear 
equations of modified Kadomtsev-Petviashvili in obtaining solutions of the bilinear 
equations is also pointed out. The different types of localized coherent structures driven 
by time-dependent boundaries are presented in section 4. In section 5, the behaviour of 
IE I in the background of stationary boundaries is analysed. Importance of boundaries in 
generating the various types of localized structures is also discussed. Finally, the results 
are summarized in section 6, where the importance of arbitrary function in expressing 
solutions is pointed out. 

2. Bilinearization of IE 

By making a stereographic projection of the spin of unit sphere on a complex plane, the 
spin components can be written in terms of the stereographic variable u [9] as 

S+ = S 1 + i S 2 = -^ T2 , Ss = ^-MJ (2.1) 
1 + \oo\ 2 1 + 

and eq. (1.1) takes the form 

2lo* 

iuo t + oo xx + o- 2 u yy - {oo 2 x + cr 2 ujy) - i<j) y 0J x - i4> x Uy = 0, (2.2a) 



1 + \U) 

Aio 2 

:i + m 2 ) 



We find that this form is more convenient for further analysis as discussed below. 



2.1. Painleve singularity structure analysis 

We can confirm the integrability nature of eq. (2.2), by performing a Painleve analysis 
of it. Denoting oo and oo* by F and G respectively, we rewrite eq. (2.2) as 

(1 + FG)[iF t + F xx + a 2 F yy - i{<P y F x + <f> x F v )\ - 2G(F 2 + a 2 F 2 ) = 0, (2.3a) 
(1 + FG)[-iG t + G xx + a 2 G yy - i(4> y G x + <f> x G v )\ - 2F{G 2 X + a 2 G\) = 0, (2.36) 
(1 + FG) 2 (<p xx - a 2 <j) yy ) = Aia 2 (F y G x - F x G y ), (2.3c) 

where eq. (2.3b) is the complex conjugate of eq. (2.3a). In order to carry out a singularity 
structure analysis of eq. (2.3), we effect the following local Laurent expansion for 
each dependent variable in the neighbourhood of a noncharacteristic singular manifold 

*l>(x,y,t) = 0: 

oo oo oo 

F = ij m Y,F ] (x,y,t)^, G = ,j> n Y,G i (x,v,t)V, = ^ p ^^(x jl /,t)^. (2.4) 

3=0 3=0 3=0 
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We now substitute eq. (2.4) into eq. (2.3) and look at the leading order behaviour of ip. 
Here we come across two different branches, one at m = 0,n = — l,p — and another 
one at m = —1, n = 0,p = 0. In both the cases, F , Go and O are found to be arbitrary- 
functions of x, y and t and the resonances are found to occur at j = —1,0,0,0,1,1. 
Further analysis confirms that two arbitrary functions occur at the resonance values of 
j = 1, without the introduction of any movable critical manifold, while the remaining 
coefficients can be expressed in terms of the earlier ones. Hence eq. (2.2) passes the 
Painleve test and confirm its integrability nature. 



2.2. Bilinearization in laboratory coordinates 

To construct a formal Backlund transformation, we truncate the Laurent series at the 
constant level term, that is (for the case m — — 1, n — and p = 0) 

F = F ^- 1 + F 1 , G = G , = 0o- (2.5) 

We can now construct the bilinear form of eq. (2.2), by considering 

F l = 0. 

Let F = g and ip — /, then oj — j. Under this transformation co = j,, where g and / 
are complex functions of x, y and t, eq. (2.2) can be written in terms of the Hirota's 
D-operators (which are defined as D % x D^ y D^a(x, y, t) ■ b(x, y, t) = (d x — d x >) l (d y — d y i) j 
(d t - d t >) k a(x, y, t)b(x', y', t') \ x=x , , y=y , tt=t ,) as 

(iD t -D 2 x -a 2 D 2 y )(f* -g) = 0, (2.6a) 

(iD t -Dl-a 2 Dl)(f*.f-g*.g)=0, (2.6b) 

D x (D x (f* -f + g*-9))- (/7 + 9*9) = v 2 D y (D y (f* • / + 9* • g)) • (/7 + 9*9), (2.6c) 

so that 

, . 2 D v(f* ■ f + 9* ■ 9) , or k 

<p x = -2ia !L 77rr - —s , 2.6cf) 

(f*f + g*g) 

Note that eq. (2.6c) arises due to the compatibility condition <p xy = <p yx and it is 
biquadratic. 



2.3. The line solitons 

Now, the construction of the line soliton solutions to IE becomes standard and we briefly 
indicate their forms. One expands the functions g and / as power series in the arbitrary 
parameter e as follows, 

oo oo 
n=0 n=l 
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Substituting these expansions into eqs. (2.6a)-(2.6c) and equating the coefficients of 
various powers of e, we get the respective following system of equations from (2.6a), 
(2.6b) and (2.6c): 

e 2n+1 : (id t + d? + a 2 d 2 y )g 2n+1 = - £ L>'(/* fe • g 2m+1 ), (2.8a) 

k+m=n 

e 2n : idtUL - hn) - {dl + a 2 d 2 y )(f; n + f 2n ) = 

D'( (9*2n 1+ l ■ <?2n 2+ l) " £ ' hm.) ) , (2.86) 

\ni+ri2=n— 1 mi+m2=n / 

(^-^ 2 )(/* n -/ 2n ) + J D"( £ (/ 2 * ni • / 2n2 ) + J2 (9*2^-92^)) + 

\ni+ri2=n— 1 mi+m2=n-l / 

\ D x (f*2n-2)x ~ f{2n-2)x + £ Arg^+l ' fl^na+l) ~ O" 2 -Dy(/(* 2n _2)j / ~ f{2n-2)y+ 
{ ni+ri2=n— 2 

Yl D y9*2n 1+ 1 ■ 92n 2+ l) > ' (/ 2n -2 - / 2 n-2 + £ 5 2m +l ' #2n 2 +l) = (2.8c) 
ni+ri2=n—2 ) n±+n2=n—2 

with £>' = iA — D 2 X — a 2 D 2 y , D" = D 2 X - a 2 D 2 y and /„ = 0. 

2.3.1. 1-soliton solution In order to construct the exact iV-soliton solutions (N-SS) of 
eq.(l.l), we make the ansatz 

N 

9i = Yl exp *i ' Xj = l i x + m iV + n i t ^ ( 2 - 9 ) 
j'=i 

where /j, and are complex constants. As an example, we write the forms of g and 
f 2 with the help of eq. (2.8) for N = 1 as 

£i = Mexpxi, xi = /ix + miy + nit, / 2 = exp 2(xlr + ip), (2.10) 

where 

2 2 /2 

Xi = XiR + iXi/> n i = + cr 2 "^ 2 ) and exp 2ip = ai+ ,* ) 2_"2 1 (m 1 1+m » )2 MM* and M is an 
arbitrary complex constant. Now we distinguish the two cases a 2 = +1 and a 2 = — 1 
Case(i): a 2 = +1 (IE I) 
With the choice 

M (li + ll) + i(mi + ml) 
ml — ill 

the corresponding 1-SS of IE I takes the form 

s + = 2E ^R miR + m iR l11 + L ) exp ixusechxii? (2 11a) 

A + 2B tanh xir + C tanh 2 xlr 

5 3 = 1 ^ " ^ W ^ 2 , (2.116) 

A + 2 J BtanhxiR + Ctanh 2 xii? 

where 

-E = (h R + im 1R )(l 2 1R - m 2 1R ), 
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L = i(m u ml R - lf R ) + (m\ R + hil\ R + i{mul\ R - l 1R m 2 1R )) tanh^iR, 

A = 21 ir ~ 2 ^Ir(^r + IIr™u) + ™\rI\i + 2l iR l nm 3 lR + m 2 u m\ R + Zl\ R m\ R - m 6 1R , 

B = l m m m( m2 iR + l li + m li + 1 \r) + ImhimiR + ki(m\ R - l\ R ) - h R m u m 4 1R 

C = 2m\ R + 2l u ll R m\ R + l\ R l 2 u + m 2 u l\ R + Zm\ R l\ R - 2l\ R m\ R - 

2m u l\ R m\ R - l\ R . 

Case(ii): a 2 = -1 (IE II) 
Choosing 

M= (h + ll) + t(mi + ml) 
ll — %m\ 

the 1-SS of IE II takes the form 

s+ _ 2 (Iir + im 1R ) + (m u - ilu + (l 1R + im 1R ) tanh Xir) exp ixusech-XiR ^ 
1 ir + ^1/ + 2(liRm u - lum 1R ) tanh xir + (m 2 ^ + m 2 ^) tanh 2 Xir 

g 3 — I 2 ^R ± m\ R )sec\i 2 xiR ^ ^ 

1 ir + l ii + 2(/ii?mi/ - lum 1R ) tanh xlr + (raf^ + mf 7 ) tanh 2 xlr 

2-soliton solution To construct 2-SS, we take iV = 2 in eq. (2.9). Then, g 1 takes 
the form 

9i = expxi + expx2- (2.13) 
Substituting (2.13) in (2.8) and after some calculations we obtain 

f 2 = M n exp (xi + xl) + M 12 exp ( X 2 + xl) + M 21 exp ( X i + X* 2 ) + 

M 22 exp( X2 + X ;), (2.14a) 
g 3 = L 1U exp (xi + Xi + X2) + L 122 exp (xi + X2 + X2), (2.146) 
/ 4 = Kexp(xi + Xi + X2 + x;), (2.14c) 
where 

^ 2 ~ ^ 
+ U 2 - o- 2 (m* r + m s ) 2 



M r S = ; —0 huzh i ~ x 2 , r,s = l,2, 



_ (^-/^((/.-/^-^K-m,) 2 ) . 9 

rst ((Z r + Z s *) 2 - a 2 (m r + m*) 2 )((/ t + I*) 2 - a 2 (m t + m|) 2 ) ' " " """ 



K _ ((11 ~ & - * 2 K ~ ml) 2 ){{h - l 2 f - u 2 [m x - m 2 ) 2 )P 
((h + llf - u\m 1 + m\f)({l 2 + llf - a\m 2 + m\) 2 ) 

P 



a 2 m\ — l\ a 2 ml — l\ 



(h + / 2 *) 2 - d 2 (mi + m* 2 ) 2 (l 2 + / 2 *) 2 - a 2 (m 2 + m*) 2 

Making use of the fact that now u = j = 1 + 9 J^+ 9 Jf 4 ; an d the relations (2.1) for the spin 
variables in terms of u>, the spin two soliton solution can be written explicitly. 
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2.3.3. N-soliton solution Finally, by taking gi as in eq. (2.9) and extending the above 
procedure, one can obtain the iV-SS as 

2N 2N 

g= expiY&ihtiViVj + ^VibCi + ijii)]}, (2.15a) 

fij=0,l i<j i=l 



where 



and 



where 



log(<r 2 m 2 -11) for i = N + 1, . . . , 2N, 
otherwise, 



N 



N 



means 
m=o,i i=i 



^2^i= l + Y^i+N, 



i=l 



' 2N 2N 

f = Yl ex p{53^( i 'iWj+5^/ i <[xi+V' / (*)]}, 

Hj=0,l i<j 



i=l 



(2.156) 



Xi = Ux + vrciy + + x»(0), 

n . = z(/ 2 + a 2 m 2 ), Xi+iv = Xi > W = l*, m i+N = m*, n i+N = n* 



<t>{hj) 



- \og((k + Ij) 2 - a 2 (rrii + rrij) 2 ) for i = 1, . . . , TV and j = TV + 1, . . . , 2N, 
log((/j — lj) 2 — a 2 (mi — rrij) 2 ) for i — 1, . . . , N and j — 1, . . . , N, 



log(cr 2 m 2 - / 2 ) fori = l,...,JV, 







otherwise, 



' N N 

Y means ^ to = ^2 IJ ' i+N - 

fii=0,l i=l i=l 



3. Bilinearization and solution of IE I in terms of light cone coordinates: 
The role of boundaries 

In the case of most of the (2+1) dimensional integrable nonlinear evolution equations like 
Davey-Stewartson I, modified Kadomtsev-Petviashvili, and Nizhnik-Veselov-Novikov 
equations [10-13], the bilinearized forms are first transformed into systems of linear 
partial differential equations(pdes) while using the Hirota's direct method. By solving 
these linear pdes, one can construct the line solitons(as we have done in section 2) 
which are localized everywhere except along particular lines. Looking at the nature 
of line solitons, the presence of two nonparallel ghost solitons [10] (solitons which are 
visible only in the absence of the physical field) are identified. As a dromion is the two 
soliton solution made out of two nonparallel ghost solitons, they can be embedded in 
the two soliton solution to generate a (1,1) dromion. It implies that the solution consists 
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of one bound state each in the x and y directions. This can be directly extended to 
generate multidromions. If we follow the same procedure for IE, we see from eq. (2.11) 
or eq. (2.12) that the two nonparallel ghost solitons are absent here. Hence we cannot 
construct dromions from the bilinearized form (2.6). 

In the following, we analyse the IE I in a different frame of reference consisting of 
light cone coordinates £ and rj, which are defined as 

£ = \(v + v), v = \{y-x). (3.1) 

Correspondingly, eq. (1.1) takes the form (after rescaling — | — ► t';t' — ► t) 

S t (Z, r}, t) = S A (% + S m ) + feSz - (f) v S v , (3.2a) 
(j) (v = S ■ S ( A S v . (3.2b) 
The form of eq. (3.2b) suggests that one can redefine the scalar field variable <fi as 

0(£, 77, t) = $(£, 77, t) + J mi (Z, t)d£ + J m 2 (rj, t)drj, (3.3) 

where mi((, t) and 1712(1], t) are the boundaries (arbitrary functions in the indicated 
variables). Eq. (3.2) can therefore be rewritten as 

§t(£, 77, *) = S A (% + S m ) + ($ e + 77H (£, t))§z - (% + m 2 ( V , t))S v , (3.4a) 
= S ■ A §„. (3.46) 

In terms of the stereographic variable oj given by eq. (2.1), eq. (3.4) takes the form 
2uj*(uj} + uj 2 ) 

iu) t + uzz + u m - + mi (f,*))<^ - «($r, + m 2 (r],t))uj v = 0, (3.5a) 

4i 

$ ^ = (i + 1^12)2 - ^K)- ( 3 - 5& ) 

After introducing the transformation u = j, the bilinear representations of eq. (3.4) can 
be written as 

(iD t - D\ - D\ - %m x (t t)Dz + im 2 (rj, t)D v )(f* ■ g) = 0, (3.6a) 
(iA " D\ -D\- im x (t, t)Df. + im 2 (^ t)D v )(f* -f-g*-g)=0, (3.6b) 

Ds(D v (f* .f + g*.g)). (f*f + ^) = A(A(/* • / + 9* • 9)) ■ (f*f + 9*9): (3.6c) 

so that 

$£ = ■ ^/;« , ;») , ,3.6.) 

(f*f + 9*g) 

$ ^ = _2z — I * \ — • 3 - 6e 
(/*/ + 0*0) 

One observes now the explicit introduction of the boundaries mi(£, £) and m 2 (r],t) into 
the bilinearized form (3.6), which turns out to be crucial to obtain localized solutions. 
Also eq. (3.6c) is a consequence of the compatibility between eqs. (3.6d) and (3.6e). 
It may also be noted that a similar introduction of arbitrary functions in the bilinear 
form was necessiated for the (2+1) dimensional sine-Gordan equation in order to obtain 
localized solutions [13]. 
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Expanding now the functions g and / as power series, 

9 = egi + e 3 #3 + e 5 #5 + • • • , (3.7a) 
f=l + e 2 f 2 + e 4 / 4 + e 6 /e + • • • (3.76) 

and substituting them into eqs. (3.6a)-(3.6c), we will obtain the following set of linear 
pdes by equating the various powers of e: 

e : (iD t - D\ - D\ - im^, t)D ( + im 2 (r), *)£>„) (1 ■ 9l ) = 0, (3.8a) 
e 2 : (iD t - D\ - D\ - im^, t)D^ + tm 2 ( V , t)D v )(l ■ f 2 + / 2 * ■ 1 - g{ ■ 9l ) = 0, (3.86) 
Ds(D v (l ■ f 2 + /* • 1 - g{ ■ 9l ) ■ 1) + D v (Ds(l ■ f 2 + /* ■ 1 - g\ ■ 9l ) • 1 = (3.8c) 

and so on. Let us solve the above linear equations to obtain the solutions to IE I. 
Eq. (3.8a) can be rewritten as 

Wit + 0i« + girm ~ imi(€, t)gi£ + im 2 (rj : t)g lr) = 0. (3.9) 
Let us try a variable separation, by postulating 

9i=p(S,t)q(v,t), (3.io) 

where p and q are complex functions of the indicated arguments. Now, eq. (3.9) becomes 

q(ip t +Ptt - im^tjpt) +p(iq t + q m + im 2 (rj : t)q v ) = 0. (3.11) 

The above equation suggests that we should have 

iPt - imi(£,t)pt = hp, (3.12a) 
iqt + q m + im 2 (rj, t)q v = - faq, (3.126) 

where k\ is a constant. Redefining p = pexp —ikit and q = qexpikit and dropping the 
hats, eq. (3.12) becomes 

m+Pti ~ imi(Z,t)pz = 0, (3.13a) 
iqt + Qrm + i^iv, *)?»7 = °- (3.136) 
Simplifying now eq. (3.8c), we will get 

f'kn ~ far, = 9*i9Uv - 9i9^r,- (3-14) 

Let us write f 2 = f 2 R + if 2 i, where f 2 R and f 2 i are the real and imaginary parts of the 
function f 2 , respectively. Eq. (3.14) now becomes 

fcmv = ^(9i9iz v ~ 9i9Hv)- 
On carrying out integrations with respect to £ and 77, we get 

hi= kil S {919 ^ ~ 9 * 9i ^ )didr} + k ^ i) + r ^ t} ' (3 - i5) 

where h(£, t) and r(rj, t) are arbitrary functions in the indicated arguments. Substitution 
of g\ and f 2 in eq.(3.8b) gives rise to 

V 2 / 2 « = f2m + f-zRm = hit - mi(£, t)f 2 i£ + m 2 (r], t)f 2Ir) - {pg>\qq* + PP*q v q*), (3.16) 
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which is nothing but the Poisson's equation in two dimensions for f 2 R. For the given 
boundaries mi(£, t) and m 2 (r), i), once the eqs. (3.13a) and (3.13b) are solved for p and q 
respectively, the right hand sides of eqs. (3.15) and (3.16) can be found out. Hence they 
can be solved to obtain f 2I and f 2 R- At this point, we may point out that exactly the 
same kind of linear equations (3.12a) and (3.12b) have appeared both in the case of 1ST 
method [7] and binary Darboux transformation method [8]. Now these equations have 
appeared from a different perspective, namely from the point of view of bilinearization 
method. 

Now the general solution for oj can be written as 
g p(£,t)q(r},t) 

J 1 + J2R + IJ2I 

where the functions p, q, f 2 j and f 2 R are to be determined from equations (3.13a), 
(3.13b), (3.15) and (3.16), respectively. Here the main task lies in solving eq. (3.13) 
which is also associated with the linear problem of the following modified Kadomtsev- 
Petviashvili equation(mKP) [14]: 

Ut + Uxxx + 6u 2 u x — 12d~ l u yy + Ylu x d~ x u y = 0. (3.18) 

The mKP eq. (3.15) is equivalent to the compatibility condition for the linear system, 

(2id y + d 2 x + 2iud x )ip = 0, (3.19a) 

(id t + Adl + I2iud 2 x + (6iu x + V2d~ x u y - 6u 2 )d x + a)ip = 0, (3.196) 

where a is an arbitrary constant. By comparing eq. (3.13) and eq. (3.19), it is seen that 
the proper solutions of eq. (3.13) are obtained by dropping the time dependence in the 
various types of solutions such as line solitons, line lumps and line breathers of eq. (3.19) 
and changing y — > t. Since we have two independent problems in eq. (3.13), large 
number of classes of exact solutions of IE I, which include rationally localized, rationally- 
exponentially localized and fully exponentially localized structures, are possible. We 
present the details of some of them in the next section. 

4. Localized coherent structures of IE I: Time dependent boundaries 

4-1. Lump-lump boundaries: Algebraically decaying structures 

In this section, we choose specific forms of the boundaries rai(£, t) and m 2 (r),t) in 
eq. (3.8) to obtain localized solutions. For the algebraically decaying boundaries of the 
form (for easy comparison, we follow the notations of ref. [7]) 

mi(£,t) = 2a 5-, (4.1a) 

2d 

m 2 ( V ,t) = - (4.16) 

(V'j + c 2 ) 2 + £- 

where a, j3, C\ and c 2 are real constants, the solutions to eqs. (3.13) are obtained as 
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expi(| — -L) 

' p ~ 2 2 

Substitution of eq. (4.2) in (3.15) leads to 

1 flg-f + + + <*) 

2 «^((e-f + c 1 ) 2 + f)(( ?? -f + c 2 )2 + ^)' 

wherein t) = and r(rf, i) = and substituting the values of p, q and / 2 j in eq. (3.16) 
and solving the resultant Poisson's equation for fan, we get 

, _ 1 (^f +C!) + fo-f +C 2 )-f 

+ C i) 2 + f)((^-f + c 2 ) 2 + ?)' 

Therefore, 

A = & + lA ' = a/3K -t + Cl -f) ( ,-| + C2 -f)- (45) 

Hence, from eq. (3.17), we have 

'1 L _l 1 L 

w = i ~~77 o,. . 7777 777" 777- ( 4 - 6 ) 



eX PHa a 2 + g g 2 - 



2 / 

Using eq. (2.1), the spin components corresponding to this form of a; are given below: 
c+ _ (l4 + (e-f + C!-f)(r/-f + c 2 -¥))expz(|-^ + |-^) 
^ ~~ z A 2 + B 2 ' 



where 



and 



*3 = 1 - -^TTT^^ (4-76) 



ap ct p 4 



The solution (4.7) decays as j- as £ 2 + r/ 2 — ► oo and moves with the velocity 
V = (Ye, V v ) = (|, |). A snap shot of S 3 is plotted in fig. (1) 

4-2. Soliton-soliton boundaries: Exponentially decaying dromions 

Let us now choose the boundaries mi(£, t) and m 2 (r/,t) as the following line 
solitons(again using the notations of [7]): 

Ml exp Ielp 

mi(£,t) = lM|2 9 J mP 7 -^, (/i = + W (4.8a) 

(l + 5fexp^) 2 + ex P ^' 

Ml eX p 2Al^ 
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Figure 1. Algebraically decaying lump structure 



where £ = £— j^f t + £oi ?7 — ? ? — j^f^ + ^o an d /Ar, A*/) -\r, A/, £ and i]q are real constants. 
The suffices R and / are used to denote the real and imaginary parts respectively. On 
solving eqs. (3.13), (3.15) and (3.16), we get the following expressions for p(£,t), q(r},t) 
and f(C,rj,t): 

expi(£ — y 
( 1 + ?7 ex P #0 

*(v,t)= . e ^-il y (4.96) 
(1 + ^expp^r/) 

_ i a + if ex p fpj) ex p wl + (j + g ex p ) ex p , 4jn 

4 ll + ^exp^ |l + Aexp^ 

f a + g -p g + g -p fgg - -p 2 m + (49d) 

4|l + ^exp^ |l + Aexp^ 

Consequently, we have 

exp ,(i + * ) 

cu = V „ " * ^ . (4.10) 

( 1 + & ex P^)( 1 + S ex P^) + i 

The spin components corresponding to this form of uj are given below: 
s+ _ W + g ^P ^0(1 + g exp + |) exp i(l - 4, + j - & 
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Ss = l IH 2 1 " , (4.116) 

where 

= 1 + 7 + — exp — £ + — exp — n + {\ R fi R - Xi^i) exp(^£ + — ^rj), 

4 m \/j,\ 2 A/ \X\ 2 Xnn \ii\ 2 \x\ 2 

D = exp + exp —n + (— + — exp — £ + — r? . 

I a* I 1^1 Ar A 4 / I a* I 1^1 

The solution (4.11) decays exponentially (see fig. (2)) in all directions on the plane £, 77 
and moves with the velocity V = (V 6 , V v ) = (gf , Qf). 




^.5. Lump-line soliton boundaries: Rationally- exponentially decaying nature 

Taking the choice of the rational lump as the boundary rai(£, £) and of the line soliton 
as the boundary m 2 (r?,£), that is 

2a 

m 1 (£,t) = ^ j", (4.12a) 

(e-f + ci) 2 + f 

eX p 

and proceeding as before, we can find the solution of u> as 

w = n ex^(i-^ + f-^) 

£ + + Cl -f)(l + gexp^) 
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The spin components in this case are 

s+ = 2(^ + (e-f + C i + f) + (l + iexp^))exp^-^ + f-^) 

£2 + F 2 + exp 2A^ 



2exp^?) 

S 3 = l |A| ' , (4.145) 

E2 + F 2 + exp ^ V ' 

where 

# = (? + ci) - - exp — 77 + — (f + ci) exp —77 + — 

a 2 \X\ Z A/ a \X\ Z 2a 

a X R 2Aj ^ , 2rj 2A/ A 

F= 2 (1+ A7 eXP ^ ) + ( ^a +Cl)eXP ^ 

The above solution (4.14) decays rationally in the £ direction and exponentially in the 
77 direction on the plane £, 77 and it moves with the velocity V = (V^, V v ) = (-, ?Af ) (see 
(3))- 




Figure 3. Rationally-cxponcntially deacying structure 

We can take the other choice also, that is mi(£, £) as the plane soliton boundary 
and 777,2(77, t) as the rational boundary. Replacing £ by 77, a by /3, 77 by | and A by /i in 
(4.12) and (4.13), we will get the solution corresponding to this case. 



5. Stationary boundaries 

In the case of stationary boundaries, that is 777,1 (£, t) = mi (£) and 771,2(77, = m 2(v)i the 
function can be written as 

<& = pWfa)T(f) (5.1) 
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where p', q' and T are complex functions of their indicated arguments. Substitution of 
this expression for g 1 in eq. (3.9) leads to 

ip'q'T t + q'T(p^ - im^)^) + p'T{q' m + im 2 (r])q' v ) = 0. (5.2) 

Eq. (5.2) suggests that 

p'x - imi(£)z4 = ~ X2 P' ( 5 - 3a ) 
q' m + irn 2 (v)q' v ) = -A'V (5.36) 
T t + i{\ 2 + X' 2 )T = 0, (5.3c) 

where A and A' are some complex arbitrary parameters. Eq. (5.3c) can be solved 
immediately and it gives 

T = exp-i(A 2 + A' 2 )t. (5.4) 

Rescaling = p(0 ex P t / m i(£)^£ an( i tfiv) = <j(v) ex P "y / ^{jfjdn and then 

dropping the hats, eqs. (5.3a) and (5.3b) take the form 

pU + ^ + ^ + ^^f)p' = 0, (5.5a) 

4 + (A ,2 + x- 2! T^' = - (5 - 5&) 

The forms of f 2 i and /2R remain the same as in eqs. (3.15) and (3.16), respectively, 
except that here we have to replace p by p' and q by q' and the solution to uj is also 
given by (3.17). Again the above linear equations are exactly the same as the ones 
obtained in ref. [7] through 1ST formalism, which are however now obtained from the 
bilinear form of IE I. 

It is seen that the problems (5.3a) and (5.3b) are gauge equivalent to the spectral 
problems of the Schrodinger type 

<j>xx + (A 2 + u 2 (x) ± iu x (x))(p(x, A) = (5.6) 

with the very special potential [14] 

V(x) = -(u 2 (x)±iu x (x)), (5.7) 

where u(x) is a real valued function. In this case also, the class of coherent structures 
admitted by (5.6) and hence by IE I are very rich. 

5.1. Constant boundaries: Lump solution 

If the boundaries are constant, that is mi(£, t) = mi and m 2 {r\, t) = m 2 , then going 
through the various steps indicated above, the function u is found out to be 

expi((m 1 + !)(£-!) + (-m 2 + I) (77 - f )) 
u = j ^ — — . (5.8) 
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The corresponding spin components are 

+ (-2 + + jSA) exp i(K + + ("^ + J) - |)) 
5+ = - f r , (5.9a) 



(l - f -, ttt^ + r-^) ^f ) 2 + 1 



5 3 = 1 j s — — • (5.96) 

(! - ((^TX)2 + (^rp J ^) 2 + 1 

In this case, the solution decays algebraically as ( - V 1 ) -1 as £ 2 + rf — > oo. 



5.2. Boundaries are absent: Line solitons 

If t) = m 2 (rj,t) = 0, then on solving equations (3.13), (3.15)-(3.17), we will get 

the following results: 

p= Aexp(l£ + il 2 t), (5.10a) 
q = B exp (mr/ + im 2 t), (5.106) 

hi = k 1 e X p2(l R £ + m R ri-2(l R l I + m R m I )); h = —AA*BB* ( ^Jre^ ) ' (5 ' 10c) 
f 2R = k 2 ex P 2(l R ^ + m R r ) -2(l R l I + m R m I )); k 2 = AA*BB* ( ^j"^*™* ) , ( 5 - 10c O 

ABexp(/£ + m?? + i(/ 2 + m 2 )f) lm AA*r>r>* /r in x 

a; = ; L — ; rr-; fc 3 = ; ijii , (o.lOe) 

1 + k 3 exp 2(/ jR £ + m R r) - 2{l R l I + m R m/)) 4/^mR 

which is nothing but the line soliton solution similar to the line solitons obtained in 
section 2 except for the difference in constant values. Here A and B are complex 
constants. So, if the boundaries are absent to start with in the IE I, eq. (3.2) or its 
bilinearized version (3.6), we obtain line soliton solutions only. From this, one can 
appreciate the important role of boundaries to generate localized coherent structures in 
the bilinearized version of IE I. 



6. Conclusions and Discussions 

If we look at the eqs. (3.13), (3.15), (3.16), it is seen that we are having four equations 
for six unknowns viz., p(^,t),q(rj,t),m 1 (^t),m 2 (r],t), / 2j r(C, ??, and f 2 i(£,r], t )- Hence 
any two of them can be treated as arbitrary functions. In this case, one can manipulate 
these arbitrary functions to get a large number of solutions similar to the solutions 
generated by Lou and his coworkers for equations such as Davey-Stewartson model, 
Nizhnik-Novikov-Veselov system and dispersive long wave equation [15]. For example, 
if we choose t) = exp xi and q(rj, t) = exp X2, where Xi — lx+uit and X2 = my+uj 2 t, 
we will get the line solitons. It should be noted that we can interpret our analysis on 
localized coherent structures done in section 4 entirely in a different way. That is, choose 
the values of t), and q(rj, t) given in eqs. (4.2a) and (4.2b), then solve the eqs. (3.13a), 
(3.13b), (3.15) and (3.16) to obtain the values of mi(£, t), m 2 (r), t), f 2 i(£,, rj, t) and 
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/2i?(Ci Vj t) respectively. Once again we arrive at the same algebraically decaying 
structure of IE I given in eq. (4.7). Similar interpretation holds good for generating 
localized structures of other types also. 

On the otherhand, if we take the more general solution of mKP I equation and 
by solving eqs. (3. 13), (3. 15) and (3.16), we can generate the multi-soliton (multilump/ 
multidromion) solutions for the IE I. For example, we take 

mifo t) = — ^ , (6.1) 

where Xj = £ — — + S i: i — 1,2, which describes the scattering of two lumps of the form 
given in eq. (4.1a) and m 2 (r/, t) given in eq. (4.1b). On solving eqs. (3.13), (3.15) and 
(3.16), we can generate a (2,1) lump solution for IE I. Hence by appropriately choosing 
the boundaries rai(£, t) and 7712(77, t), one can generate multidromions or multilumps 
and in general multisoliton structures. 

To conclude, in the case of IE also the Hirota method has proved to be a 
straightforward but powerful tool to obtain different kinds of localized and other 
solutions. This has become feasible by appropriate bilinearization procedure in the 
light cone coordinates. We have also reported how the explicit forms of various localized 
solutions can be deduced. 
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